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Abstract: 

Throughout  the  comprehensive  research  in  the  past  three  years,  the  unique  properties 
of  Type  II  InAs/GaSb  heterojunctions  were  utilized  for  the  realization  of  novel  uncooled 
infrared  photodetectors  with  higher  operating  temperature,  detectivity  and  uniformity 
than  the  commonly  available  infrared  detectors.  We  have  successfully  demonstrated 
uncooled  photon  detectors  in  the  long  wavelength  infrared  range  from  8  pm  to  12  pm 
based  on  Type  II  InAs/GaSb  superlattices.  Photoconductors  were  grown  on  GaAs 
substrates  or  on  GaSb  substrates  with  AlSb  insulating  layer.  Photodiodes  were  grown  on 
GaAs,  GaSb  substrates.  Wet  chemical  etching  and  dry  etching  techniques  have  been 
developed  successfully.  For  photoconductors  at  room  temperature,  a  detectivity  of 
1 .08x1 08  cmHz1/2/W  at  1 1  pm  and  an  effective  carrier  lifetime  of  26  ns  have  been 
achieved.  For  photodiodes  at  room  temperature,  we  have  obtained  a  detectivity  of 
1.2xl08  cmHz1/2/W  at  7  pm,  and  R0A=1.36xlO'2  Gem2,  better  than  HgCdTe  (~10‘4  Gem2 
).  Detector-laser  measurements  have  shown  a  response  time  of  less  than  68  ns  for  these 
detectors.  However,  some  properties  are  still  far  from  the  focal  plane  array  application 
requirements.  The  detector  performance  still  needs  to  be  improved  in  the  continuation  of 
the  program,  especially  the  detectivity  and  RqA. 


I. 


Introduction 


Infrared  (IR)  detectors  are  the  critical  parts  of  many  modem  medical,  industrial, 
military,  and  scientific  instruments  and  systems.  Some  examples  are  thermal  imaging 
systems  for  early  detection  of  breast  cancer,  non-invasive  glucose  monitoring,  fast  engine 
inspection  in  the  aviation  industry,  infrared  active  countermeasure,  and  laser  RADAR. 
The  demand  for  high  performance  IR  detectors  in  these  systems  has  driven  academic  and 
industrial  research  and  development  toward  exploration  of  new  material  systems  as  well 
as  device  stmctures.  Currently,  the  most  commonly  used  material  in  the  infrared  systems 
is  HgCdTe.  Although  it  is  the  most  studied  semiconductor  after  silicon  and  germanium, 
it  still  cannot  provide  the  required  performance  due  to  several  physical  limitations.  Some 
of  these  include  non-uniformity  of  the  composition  due  to  the  high  Hg  vapor  pressure, 
low  mechanical  strength  and  radiation  hardness  due  to  the  weak  atomic  bonding,  and 
high  tunneling  current  due  to  the  low  electron  effective  mass. 

Type  II  InAs/GaSb  superlattices  have  been  proposed  as  an  alternative  to  HgCdTe 
for  the  strategic  windows  of  3-5  pm  and  8-12  pm  as  well  as  longer  wavelengths.  Unlike 
HgCdTe,  Type  II  superlattices  are  constructed  from  the  III-V  material  system,  and  hence 
they  have  much  better  mechanical  properties  and  material  uniformity.  The  electron 
effective  mass  in  these  superlattices  is  higher  than  the  electron  effective  mass  in  HgCdTe 
and  InAsSb,  therefore  the  tunneling  current  is  less.  The  bandgap  of  the  superlattice  can 
be  changed  to  cover  a  wide  IR  range  (from  ~2pm  to  above  25pm)  by  changing  the 
thickness  of  the  layers  rather  than  the  material  composition.  Also,  in  comparison  to 
quantum  well  infrared  photodetectors  (QWIPs),  Type  II  superlattices  have  a  much  higher 
quantum  efficiency  while  they  do  not  require  surface  gratings  for  absorption  of  normal 
incident  light. 

Despite  such  inherent  advantages,  only  a  few  successful  IR  detectors  based  on 
Type  II  superlattices  have  been  demonstrated  recently.  The  major  difficulties  in  the 
realization  of  Type  H  detectors  are  the  design  of  the  optimum  superlattice  structure,  the 
growth  of  high  quality  material,  and  the  processing. 

The  goal  of  our  study  is  to  exploit  the  advantages  of  Type  II  superlattices  for  the 
realization  of  high  performance  infrared  detectors  at  room  temperature.  First, 
comprehensive  theoretical  modeling  has  been  developed  for  the  simulation  of  the  band 
structure  of  the  Type  II  superlattices  as  well  as  the  simulation  of  the  detectors.  Using 
these  simulation  programs,  the  optimum  superlattice  and  detector  structures  for  higher 


operating  temperatures  and  higher  quantum  efficiencies  have  been  derived.  In  parallel, 
the  optimum  conditions  for  the  growth  of  these  superlattices  have  been  determined  based 
on  their  structural,  optical,  and  electrical  characterization.  Finally,  for  the  fabrication  of 
the  devices,  several  processing  techniques  such  as  wet  and  dry  etching  of  the  superlattice 
structures  have  been  developed. 


II.  Realization  of  Uncooled  Type  II  Photoconductors  in  the  Long 

Wavelength  IR  range 

II.A.  Motivation 


II.A.l.  The  Need  for  Uncooled  Photon  IR  Detectors 

Currently  available  photon  detectors  have  low  operating  temperatures,  and  hence 
require  cryogenic  coolers.  However,  in  most  of  the  applications,  these  coolers  are  not 
desirable  because  of  their  short  lifetime  and  the  added  power  consumption,  weight, 
volume,  and  costs.  Commercially  available  uncooled  IR  imaging  sensors  use 
ferroelectric  or  microbolometer  detector  arrays.  These  sensors  are  inherently  slow  and 
cannot  detect  rapid  scene  changes  needed  for  many  applications.  Some  of  the 
applications  which  require  a  fast  detector  response  time  (x  <  30  msec)  are:  ffee-space 
communication,  proximity  fuzes,  active  infrared  countermeasure  systems,  missile 
detection/situational  awareness  for  highly  maneuvering  airborne  platforms,  LIDARs, 
gated-imaging,  and  night  vision  systems.  Thus  there  is  a  need  for  the  development  of 
high-speed  uncooled  detectors  in  order  to  meet  the  requirements  of  present  and  future 
applications. 


II.A.2.  Problems  of  Currently  Available  Photon  Detectors 

Although  photon  detectors  have  gigahertz  bandwidths,  their  high  temperature 
detectivity  is  severely  degraded  due  to  several  physical  limitations.  The  existing  infrared 
photon  detectors  can  be  categorized  as  interband,  which  are  mostly  HgCdTe  and  InAsSb, 
or  intersubband  quantum  well  infrared  detectors(QWIP).  There  are  some  fundamental 
limitations,  namely  fast  Auger  recombination  rate  in  the  interband  detectors  and  high 
thermal  generation  rate  in  the  intersubband  detectors,  which  drastically  decrease  their 
performance  and  ability  for  high  operating  temperature.  Moreover,  the  difficulty  of  the 
growth,  nonuniformity  due  to  high  sensitivity  to  the  composition,  and  large  tunneling 
currents  in  HgCdTe  and  the  required  sophisticated  processing  for  normal  incidence  light 
coupling  in  n-type  QWIPs  are  the  other  drawbacks  of  the  currently  available  IR  photon 
detectors. 


H.A.3.  Advantages  of  Type  II  Superlattices  for  Uncooled  IR  Detection 

In  comparison  to  HgCdTe,  the  higher  effective  mass  of  electrons  and  holes  and  the 
slower  Auger  recombination  rate1,2  lead  to  lower  dark  current  and  higher  operating 
temperature  in  Type  II  superlattices.  Another  advantage  of  a  Type  II  superlattice  is  the 
possibility  of  bandgap  engineering.  Unlike  bulk  material  or  Type  I  superlattices,  one  can 
modify  the  energy  of  the  conduction  and  valence  minibands  of  a  Type  II  superlattice  with 
a  high  degree  of  freedom.  Recently,  II- VI  HgTe/CdTe  and  III-V  InAs/GaxIni-xSb  Type  II 
superlattices  have  shown  very  promising  results  in  the  long  wavelength  ranges  ’  . 

II.B.  Material  System 

InAs,  GaSb,  and  AlSb  and  their  related  compounds  form  a  nearly  lattice  matched 
family  of  semiconductors  known  as  the  6. 1  A  family,  since  the  lattice  constants  of  these 
material  are  about  a=6. 1  A.  Figure  1(a)  shows  their  energy  gaps  and  lattice  constants 
compared  to  the  other  major  semiconductors,  while  Figure  1(b)  shows  their  energy  band 
lineups. 


(a)  (b) 

Figure  1  (a)  The  energy  gap  and  lattice  constant  of  InAs,  GaSb,  and  AlSb  and  related 
ternaries  compared  to  the  other  major  semiconductors,  (b)  The  energy  band  lineups  of 
InAs,  GaSb,  and  AlSb. 


The  following  table  shows  some  of  the  important  properties  of  the  6.1  A  family: 


Parameters@300K  (unless 
indicated  otherwise) 

InAs 

GaSb 

AlSb  | 

Lattice  constant  (Angstrom) 

6.058 

6.094 

6.136 

Density  (g/cmJ) 

5.70 

6.61  ] 

4.26  1 

Melting  point  (K) 

1215  _1 

985 

1338 

- - - jr - 

Expansion  coefficient 

4.52x10’° 

6.2x1 0’° 

4x10'° 

Thermal  conductivity(WK"'cm‘') 

0.48 

0.4  1 

~0.7 

Energy  gap  (  electron  volt) 

0.356@300K 

0.414@77K 

0.7@300K 

-~-@77K 

1.63@300K 

1.666@77K 

Electron  mobility(cmz/v.s) 

20000@300K 

35000@77K 

5000@300K 
. — (g)77K 

200@300K 
. — @77K 

Hole  mobility(cmz/v.s) 

480@300K 
— @77K 

880@300K 

2400@77K 

375@300K 
- @77K 

Electron  effective  mass 

0.023mo 

0.042mo 

0.12  1 

Hole  effective  mass 

0.42  mo 

0.4  mo 

0.98  mo 

LO  phonon  energy  at  zone  center 

29.6meV 

28.8meV 

42.1meV 

TO  phonon  energy  at  zone  center 

26.9meV 

27.7meV 

39.5meV 

Refractive  index 

3.42@10pm 

3.52@3.7pm 

4.56@0.517pm 

3.84@10pm 

3.79@2|0.m 

3.92@1.55pm 

2.995@20pm 

2.08@15pm 

3.1@10p,m 

3.3@2pm 

Dielectric  constant  (static) 

15.15 

15.7 

12.04 

Table  I.  Some  of  the  important  properties  of  InAs,  GaSb,  and  AlSb. 


The  maximum  lattice  mismatch  in  this  material  system  is  about  ±0.6%,  and  hence 
many  high  quality  structures  can  be  designed  within  the  critical  thickness  of  the  layers.  It 
is  also  possible  to  balance  the  mismatch  of  the  heterostructure  or  superlattice  such  that 
the  average  lattice  constant  of  the  structure  is  equal  to  the  substrate.  Another  advantage 
of  this  material  system  is  the  wide  range  of  energy  gaps  as  well  as  band  lineups  that  are 
available.  Superlattices  based  on  InAs,  GaSb,  and  AlSb  layers  can  be  designed  with 
tunable  bandgaps  from  about  1 .6  eV  down  to  zero  (semi-metal),  while  conduction  band 
offsets  above  1 .4  eV  and  valence  band  offsets  up  to  0.6  eV  are  possible. 


II.C.  Modeling  and  Simulation  of  the  Superlattices  Using  k*p  Method 
II.C.l.  Energy  Band  Modeling 

The  modeling  of  the  superlattice  energy  bands  is  a  necessary  step  in  the  realization 
of  Type  II  infrared  detectors  since  it  provides  crucial  information  about  the  bandgap  of 
the  superlattice  and  the  electron  and  hole  wavefunctions.  Based  on  this  information,  one 
can  design  a  superlattice  structure  for  a  specific  cut-off  wavelength  with  maximized 
absorption  coefficient  and  carrier  lifetime. 

The  basis  for  such  modeling  is  the  envelope  function  approximation  in  which  the 
band  structure  of  a  periodic  or  non-periodic  heterojunction  (superlattices  or  quantum 
wells)  can  be  modeled.  In  this  approximation,  we  assume  that  the  electron  or  hole 
wavefunctions  consist  of  two  parts:  a  periodic  part  due  to  the  regular  crystal  periodicity, 
and  an  envelope  part  due  to  the  heterojunction.  The  envelope  modulates  the  periodic  part 
and  is  assumed  to  be  much  larger  than  it. 

Assuming  a  heterojunction  of  material  A  and  B,  the  electron  wavefunction,  \|/(r), 
can  be  written  as: 

4'(r)  =  If/A'B,(r)u,(r)  (1) 

/ 


where  r  is  the  position  vector  in  real  space,  f/A,B^  (r)is  the  envelope  function  in  layer  A 
and  B,  and  u/  (r)  is  the  periodic  part  of  the  Bloch  function,  and  /  runs  over  as  many  bands 
as  are  included  in  the  analysis.  The  envelope  functions  can  be  decomposed  into 
components  that  are  in-plane  and  perpendicular  to  the  A-B  junction: 


f  (A,B) 


(r-z)=7T 


(»krr±)^  (A>B) 


(z) 


(2) 


where  rx  and  z  are  the  perpendicular  and  parallel  to  the  growth  direction  vectors,  kj_=(kx 
,ky)  is  the  perpendicular  wavevector,  S  is  the  area  of  the  sample,  and  x/A,B)( z)  is  the 
parallel  envelope  function  for  band  /.  The  main  goal  is  the  calculation  of  the  %  functions. 
The  Hamiltonian  for  the  heterojunction  is: 

H  =  ^-  +  VA(r)YA+VB(r)YB  (3) 

2m  o 

where  p  is  the  momentum,  mo  is  electron  mass,  VA(r)  and  VB(r)  are  the  atomic  potential 
in  layers  A  and  B,  and  YA=1  in  layer  A  and  zero  otherwise,  and  YB=1  in  layer  B  and  zero 
otherwise.  We  have: 

Hu/0(r)  =  (E/0(A)Ya  +  E,  0(B)YB)u/0(r) 


(4) 


where  u/,o(r)  is  the  periodic  wavefimction  and  E/,o  is  the  energy  of  band  /  at  k=0.  Now 
applying  the  Hamiltonian  to  v|/(r)  and  simplifying  it,  we  find  that  the  %  functions  should 
fulfill  the  following  for  energy  E: 

dx  =  ex  (5) 


where  X=(X/>  X™ _ )  is  an  lxN  vector  of  x/  for  N  different  bands  and  D  is  a  NxN  matrix 

with  elements  D/  „,  as  following: 


D;,m  - 


E;,0(A)Ya+E;,oWYb  + 


(B)a 


r 


h2  kx2 


2mn  2mn  dz 2 


>l,m 


tikL 

2m„ 


•</|px  |»7>- 


(6) 

2m0  dz 


where  —  and  —  are  the  first  and  second  order  derivative  operators,  l  and  m  are  two 

dz  dz 2 

different  bands,  and  5/,m  isthe  dirac-delta  function  (zero  for  h±m  and  one  for  l=m).  Now  it 
is  clear  that  matrix  D  is  just  the  k.p  matrix  of  the  bulk  materials  A  and  B,  except  that:  kz 

is  replaced  with  —and  kz2  with  ~  and  E/,0  depends  on  whether  one  is  in  layer  A  or  B. 

y  dz  dz2 

Since  the  conduction,  heavy-hole,  light-hole  and  spin  orbit  bands  have  considerable 
interaction  in  narrow-gap  Type  II  superlattices,  the  eight-band  k.p  matrix  (N=8)  was 
chosen.  The  8x8  matrix  at  ki=0  is  decomposed  to  two  4x4  matrices  of  D  andD  : 


0 

— > 

_ i 

-  p 

1 

-> 

_ 1 

i 

O 

- 1 

— > 

Q 

i 

— > 

_ i 

—  IE 

LxN 

Since  the  spin-up  and  down  are  identical  at  kj_-0  (D  -D^— D  and  %  %  X)>  eigen 

;np,  I1 


are  twice  degenerate. 

Matrix  D  is: 

Ec(z)  +  ^—  PzFPz 

0 

-f3] 

2m0 

0  Eini (z) ' 

-t^“Pz(Yi  -2y2)p2 
2m0 

0 

-|nP. 

0 

Elh  (z)  "  ITT"  P 
2m0 

|"p. 

0 

V2 

- p2 

m0 

-npz 

0 


-PzTzPz 


Eso(z)-- — PzTiP 
2m0 


(8) 


where  Id,  yi,  y2,  and  F,  are  semiconductor  parameters  and  can  be  found  in  semiconductor 
data  books  such  as  Landolt-Bomstein.  Ec(z),  Ehh(z),  Elh(z),  and  Eso(z)  are  the  energy  of 
the  conduction,  heavy-hole,  light -hole,  and  spin-orbit  split-off  bands  at  k— 0  in  the  bulk 
semiconductors.  These  are  functions  of  position  since  at  different  z  values  different 
layers  of  semiconductor  with  different  band  energies  exist.  Although  the  energy  of 


different  bands  at  k=0  is  readily  known  for  most  of  the  binary  and  ternary 
semiconductors,  the  band  lineup  of  two  semiconductors  at  their  heterojunction  needs  to 
be  calculated. 

II.C.2.  Band  Alignment  Modeling 

The  theoretical  calculation  of  the  band  lineups  at  semiconductor  heterojunctions 
has  been  a  difficult  task,  especially  since  experiments  show  a  wide  range  of  measured 
values.  Although  many  different  models  have  been  suggested5  for  such  calculations,  they 
are  not  convenient  for  our  modeling,  since  they  require  huge  computations.  Model-solid 
theory6  provides  a  simple  yet  accurate  method  for  the  calculation  of  the  band  lineups.  It 
also  does  not  require  a  posteriori  for  the  strain  effect,  since  strain  will  directly  appear  in 
the  deformation  potentials. 

Assuming  that  as  is  the  lattice  constant  of  the  substrate  and  ae  is  the  lattice 
constant  of  the  epi-layer,  the  strain  field  parallel  to  the  junction  is: 


The  strain  in  the  perpendicular  direction  is: 
8i  =  -Den 


and  the  value  of  D  for  different  crystal  orientation  is: 

pOOl  _  2  C|2  ■  pUO  _  C11  +^C12  ~2c44  ,  j-yn  _  o  cl  1  +  ^C12  ~2c44 

cn  ’  c,  |  +  c12  +  2c44  cu  +  2c12 +4c44 


(10) 

(11) 


where  cn,  C12,  and  C44  are  elastic  constants  of  the  epilayer. 

When  the  strain  is  along  [001],  the  position  of  the  heavy-hole,  light-hole,  and  spin-orbit 
bands  can  be  calculated  from: 


Ehm  —  Ev>av  +  av(2e||  +  8_l)  +2  2  ^°01 

Elh  =  EVav  +  av  (2S||  +  £1 )  ~  +  ~  §E00,  +  — 

ES0  =  Ev,av  +  av(28||  +  8x)  ^  +  J8E001 


A2o  +  A08E001  +—  (8E001)2 
Ag  +  A08E001  +  — (5E001)2 


1/2 


1/2 


(12) 


where  Ev,av  is  the  average  valence  band  energy  and  A0  is  the  spin-orbit  to  valence  band 
gap.  av  is  the  hydrostatic  deformation  potential  of  the  valence  band  and  8Eooi=2b(sj.-e||) 
where  b  is  the  shear  deformation  potential. 

The  conduction  band  can  be  calculated  from: 


(13) 


=  EVav  +  (av  +  ac)(2e||  +  ex)  + 


—  +  E 
3  8 


where  ac  is  the  hydrostatic  deformation  potential  of  the  conduction  band  and  Eg  is  the 
bandgap  of  the  semiconductor. 

II.C.3.  Energy  Band  Numerical  Simulation 

£L  Periodic  Structures 

The  Envelope  Function  Approximation  was  used  to  numerically  simulate  the  band 
structure  of  the  superlattices.  The  system  of  differential  equations  that  governs  the 
wavefunctions  can  be  constructed  from  the  relative  values  of  the  band  edge  energies  Ec, 
Elh,  EHh,  ESo,  and  the  equation  Dy=Ey  from  the  previous  sections.  This  system  can  be 
solved  with  a  proper  set  of  boundary  conditions.  A  periodic  structure  implies  a  periodic 
wavefunction  or:  %(z+L)=x(z),  where  L  is  the  period  of  the  structure  (a  superlattice  or 
multi  quantum  well).  It  can  be  shown  that  the  energy  E,  which  satisfies  the  boundary 
conditions,  for  a  two-layer  periodic  structure  is  the  solution  of: 


cos(kz  (a  +  b))  =  cos(k,a)  cos(k2b) 


2 


kiP2  +  k2Fi 


k2p,  k,p2 


■  sin(k,a)sin(k2b) 


(14) 


where  a  and  b  are  the  thickness  of  layer  1  and  2  in  the  superlattice  and  kz  is  the 
superlattice  momentum  in  the  z  (growth)  direction.  The  constants  kj  and  (A;  are  defined 
for  i=l,2  (layer  1  and  2)  as: 


k=  1 3 •  (Ec  - E)((ELH  - E)(Eso - E) - 0.5 •  8E201) 
i  A2n2(2  ■  (Eso  +  8E001  -  E)  +  Elh  -  E) 

_  3  •  (0.5  •  5Eq01  -  (E,  H  ~  E)(Eso  -  E) 
n2  (2  ■  (Es0  —  8E00i  —  E)  +  ELH  —  E) 

where  Ec,  ELh,  ESo,  and  8E0oi  are  calculated  as  described  in  the  previous  section,  n  is 
the  Kane  parameter  which  is  explained  in  the  "Energy  Band  Modeling"  section.  For  the 
heavy-hole  band: 

k  =  T-'  V^rn hh(E^7^E) 

[Mi-  _mHH 

where  EHh  is  the  heavy-hole  energy  at  k=0,  and  mHH  is  the  heavy-hole  effective  mass. 


(16) 


Figure  2  shows  the  solution  of  equation  (59)  calculated  by  a  program  written  in 
Mathcad7.  The  results  of  the  calculated  optical  transition  energies  based  on  this  model 

shows  excellent  agreement  with  our  recent  experimental  results  on  p-i-n  Type  II 

8 

photodiodes  . 


Figure  2.  The  calculated  energy  band  dispersion  of  a  InAs/GaSb  (48A/48A)  superlattice. 
kz  is  the  momentum  in  the  growth  direction,  and  d=a+b  is  the  period  of  the  superlattice. 


fx  Aperiodic  Structures 

In  the  case  of  an  aperiodic  structure,  the  above  simulation  technique  cannot  be  used 
since  the  periodic  boundary  condition  is  no  longer  valid.  However,  the  differential 
equation  derived  from  D%=E%  can  still  be  solved  numerically.  The  Runge  Kutta 
Fehlberg  method  was  chosen  as  the  numerical  technique  to  solve  the  system  of 
differential  equations,  since  it  is  proven  to  provide  high  accuracy  and  speed.  Another 
advantage  of  this  technique  is  that  if  one  uses  the  5th  order,  the  error  of  the  computation 
can  be  predicted.  This  feature  is  used  in  the  program  to  adjust  the  step  size  of  the 
independent  variable  (in  this  case  position  z)  of  the  differential  equation  to  keep  the  error 
below  a  given  value. 

The  program,  written  in  C++,  is  capable  of  calculating  the  allowed  energy  states  of 
the  heterostructures  as  well  as  the  electron  and  hole  wavefunctions  associated  to  each 


energy  state.  Figure  3  shows  an  example  of  the  simulated  electron  wavefunctions  and 
energies  of  two  states  at  the  edge  of  a  Type  II  superlattice  detector.  The  magnitude  of  the 
electron  wavefunctions  shows  that  the  quantum  wells  next  to  the  GaSb  buffer  layer  have 
a  lower  chance  of  occupation  compared  to  the  rest  of  the  wells  in  the  superlattice. 


Distance  (A) 


Figure  3.  Two  calculated  electron  energy  states  and  their  wavefunctions  in  the  active 
region  of  a  Type  II  superlattice  detector.  CB,  VB,  and  SO  are  the  conduction,  valence, 
and  spin-orbit  bands  of  different  layers. 


II.D.  Modeling  and  Simulation  of  the  Superlattices  Using  Empirical  Tight-Binding 

Method  (ETBM) 

ETBM  originated  from  the  work  of  J.  C.  Slater  and  G.  F.  Koster  in  1954.9  The 
method  was  originally  called  Linear  Combinations  of  Atomic  Orbitals  (LCAO).  The 
atomic  orbitals  are  indeed  Wannier  functions10  (for  theoretical  study)  or  Lowdin 
orbitals11  (for  practical  calculations).  Wannier  functions  and  Lowdin  orbitals  are 
orthogonal  at  different  atomic  sites,  while  actual  atomic  orbitals  are  not.  The  ETBM  use 
parameter  fitting  to  simplify  numerical  calculations,  and  is  different  from  the  ab  initio 
tight  binding  method12  which  involves  calculations  of  many  difficult  two  center  integrals 
that  are  in  general  not  practical.  We  treat  those  integrals  as  parameters  to  fit  the 
experimental  band  structure  information  and  pseudopotential  results  at  certain  symmetry 
points  in  the  1st  Brillouine  zone  for  bulk  DI-V  materials,  and  then  plug  these  fitted 
parameters  into  superlattice  Hamiltonian.  We  have  reduced  the  number  of  independent 
energy  interaction  terms  using  the  symmetry  of  the  crystal.  ETBM  can  be  used  to  treat 
III-V  semiconductors  without  magnetic  structures. 

The  ETBM  is  capable  of  calculating  the  bandstructure  in  the  entire  Brillouine  zone, 
and  simulating  arbitrary  Bravais  lattice  structures  with  basis  at  atomic  level.  This  gives 
us  the  freedom  to  do  calculations  with  very  thin  layers  in  superlattices,  including  the 
interfaces  which  are  only  one  monolayer  thick. 

There  are  different  levels  of  modeling  in  ETBM.  For  bulk  materials,  according  to 
their  accuracy,  there  are  sp3s  nearest  neighbor  model  ,  sp  s  second  nearest  neighbor 
with  spin-orbit  interactions  model14,  sp3d5s  nearest  neighbor  with  spin-orbit  interactions 
model15,  etc.  For  superlattices,  there  are  spY  second  nearest  neighbor  without  spin-orbit 
interactions  model16’17,  sp3  second  nearest  neighbor  with  spin-orbit  interactions  model  , 
etc.  Among  all  these  levels,  we  have  chosen  spY  nearest  neighbor  modeling  with  spin- 
orbit  interactions.  This  level  of  modeling  have  the  accuracy  we  need  and  with  fairly  less 
amount  of  calculations  we  can  handle. 

Another  issue  we  have  taken  into  account  is  the  stress  and  strain  effects.  The  stress 
and  strain  are  naturally  embedded  into  the  Hamiltonian  using  the  well-known  d  scaling 
rule.  The  deformation  potential  calculations  used  in  kp  modeling  are  not  needed.  The 
effects  of  the  stress  or  strain  on  the  bandstructure  can  be  calculated  directly. 

Before  we  write  down  the  Bloch  wavefunction  for  the  electrons,  we  need  to  identify 
the  simplest  unit  cell  for  the  system  to  reduce  all  redundant  calculations  and  shorten  the 
calculation  time.  For  bulk  III-V  materials,  the  unit  cell  we  choose  consists  of  two 


adjacent  group  HI  and  group  V  atoms.  For  superlattice  materials,  the  choosing  of  unit 
cells  is  illustrated  in  Figure  4. 

Assuming  the  construction  Lowdin  orbitals  for  the  bulk  or  superlattices  is. 

ft?  (r-R-f.)  (>7> 


where  n  runs  through  all  the  atoms  in  a  unit  cell,  r„  are  the  coordinate  position 
vectors  of  the  atoms  in  the  cell,  a  runs  though  all  considered  types  of  Lowdin  orbitals  (s, 
p,  s*,  etc.),  and  R  represent  the  coordinate  position  vector  of  the  unit  cell  in  the  entire 
material  we  are  considering. 


One  Superlattice  Unit  Cell 

(b) 


Figure  4  (a)  One  superlattice  example  of  [(InAs)2-InSb-(GaSb)i]n  ;  (b)  The  simplest  unit 
cell  chosen  for  ETBM  calculations. 


Now  we  can  construct  the  Bloch  wavefunction  of  electrons  for  the  entire  material. 

T-  (?)  =  £  X  X  exp(^  '  +  ^  ^  8) 

a  n=l 

where  N  is  the  total  number  of  atoms  in  one  unit  cell,  A"  are  constants,  and  RSL  runs 
through  all  the  unit  cells. 

The  bandstructure  calculation  problem  now  reduced  to  an  eigenvalue  problem  for 

^(f): 

HV-k(r)  =  EV-t(r)  (19) 

where  H  is  the  Hamiltonian  for  the  material  under  investigation,  and  E  is  the  energy 

eigenvalues  that  are  dependent  on  k  values  in  the  1st  Brillouin  zone.  Using  the 
orthogonal  properties  of  the  Lowdin  orbitals,  we  can  rewrite  the  eigenvalue  problem  in 
the  following  format: 


+ ))•  <  J ^ )^“ (r  - ^  -  ?„ >  =  ^  -  (20) 

«SI  a  »=I  “ 

For  nearest  neighbor  approximation,  RSL  runs  through  only  those  nearest  cells.  For 

the  tetrahedral  structures,  the  summation  over  RSL  will  be  reduced  to  closed  torms  and  a 
practical  set  of  equations  for  numerical  calculations  can  be  obtained.  The  only  difference 
between  bulk  material  and  superlattice  is  the  number  of  atoms  in  one  calculation  unit  cell. 
Therefore  the  matrix  size  for  superlattices  is  proportional  to  the  real  space  period  in  the 
growth  direction.  Figure  5  shows  a  general  format  of  the  Hamiltonian  matrix  for  the  bulk 
material.  Figure  6  shows  a  representative  band  diagram  in  the  entire  1st  Brillouin  zone 
along  typical  symmetry  directions  for  bulk  InSb.  The  short  blue  horizontal  lines  in  the 
diagram  are  the  energy  levels  at  certain  symmetry  points  from  the  Landolt  Bomstein  New 
Series  for  semiconductor  materials.  Figure  7  shows  one  calculated  band  diagram  for  the 
superlattice  of  [(hiAs)20-InSb-(GaSb)l3]n  along  growth  direction  and  the  in  plane 
direction. 
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Figure  5  Representation  of  the  Hamiltonian  matrix  for  bulk  13I-V  tetrahedral  structures. 
The  energy  terms  remains  as  fitting  parameters  and  f,  g  are  geometric  factors  ralated  with 
the  tetrahedral  structure.  Aa,  Ac  are  the  spin-orbit  interaction  terms  for  anions  and 
cations. 


Fitted  ETBM  Band  structure  of  InSb 


Figure  6  Representative  bandstructure  for  InSb  at  liquid  nitrogen  temperature  for  the 
entire  1st  Brillouin  zone  along  typical  symetry  directions.  The  short  horizontal  blue  lines 
in  the  band  diagram  are  energy  levels  from  Landolt-Bomstein  New  Series  for 
semiconductor  materials. 
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Figure  7  Calculated  band  structure  for  superlattice:  [(InAs)2o-InSb-(GaSb)n]n  along  the 
growth  direction  and  the  inplane  direction. 

II.E.  Modeling  and  Simulation  of  the  Devices 

In  the  previous  chapters  the  modeling  and  simulation  methods  required  for 
optimizing  the  superlattice  structure  was  described.  However,  modeling  and  simulations 
at  the  device  level  are  also  necessary  tools  for  the  optimization  of  detectors.  The  key 
issue  is  maximizing  the  detectivity  of  the  device.  It  can  be  shown  that  the  detectivity  is 
independent  of  the  device  area.  Therefore,  the  only  parameter  is  the  thickness  of  the 
active  layer.  The  thicker  the  active  layer,  the  higher  the  internal  quantum  efficiency  is 
since  the  optical  absorption  will  be  higher.  However,  a  thicker  active  layer  means  a 
lower  resistivity,  and  hence  higher  noise.  The  detectivity  of  the  device  can  be  calculated 
from  the  responsivity  and  the  noise  as  indicated  by  equation  (22) 

The  responsivity  of  the  device  as  function  of  the  SL  thickness,  t  is: 

r. (t)  __  qX.t(|ie  +  frh)VeO  ~  r)  ^  _  e~at ^ 
hcl2 


(21) 


where  t  is  the  thickness  of  the  superlattice,  q  is  the  electron  charge,  X  is  the  infrared 
wavelength,  t  is  the  carrier  lifetime,  ue  and  p,h  aic  electron  and  hole  in-plane  mobilities, 
VB  is  the  voltage  bias,  r  is  the  reflectivity  of  the  sample,  h  is  the  Planck  constant,  c  is  the 
speed  of  light,  1  is  the  device  length,  and  a  is  the  absorption  coefficient  at  X. 

The  current  noise  of  the  device  is  also  a  function  of  the  device  conductance  which 
depends  on  the  buffer  layer  and  device  thickness.  Here  we  consider  the  major  noise 
components  that  are  Johnson  and  generation-recombination  noises.  Johnson  noise  of  the 
device  is: 


in,Johnson(t)  =  ^|4kTfenlie  +  Pl^h  )*  +  Af  (22) 

where  k  is  the  Boltzmann  constant,  T  is  the  temperature,  n  and  p  are  the  electron  and  hole 
concentrations,  aB  is  the  buffer  layer  conductance,  and  tB  is  the  GaSb  buffer  layer 
thickness.  The  generation-recombination  noise  can  be  calculated  from: 


ln,GR 


(0  = 


f2qVB(He+Hh)]  f  nP 
1  l2 


n  +  p  1  +  o2t  1  j 


AdtAf 


(23) 


where  co  is  the  angular  frequency  and  co2t2«1  in  this  experiment.  The  total  current  noise 
of  the  device  can  be  calculated  from: 

—  Johnson  ”^n,GR  (24) 


A  program  in  Mathcad  was  written  which  calculates  the  responsivity  and  detectivity 
of  the  device  versus  the  superlattice  thickness.  The  input  data  for  the  program  are  the 
mobility  and  concentration  of  electrons  and  holes  (extracted  from  Hall  measurement),  the 
absorption  coefficient  (extracted  from  FTIR  measurement),  and  carrier  lifetime  (extracted 
from  noise  and  responsivity  measurements).  Figure  8  shows  the  calculated19  detectivity 
of  a  device  versus  superlattice  thickness  for  three  different  GaSb  buffer  layer  thickness. 


Thickness  of  the  SL  (iim) 


Figure  8.  The  calculated  detectivity  of  a  device  versus  the  thickness  of  the  superlattice 
with  different  GaSb  buffer  layer  thickness,  n  and  p  are  electron  and  hole  concentrations, 
pn  and  pp  are  their  mobilities,  a  is  the  optical  absorption  coefficient  and  x  is  the  carrier 
lifetime. 


II.F.  Growth 

After  the  design  of  the  superlattice  layer  thickness  for  the  target  wavelength  and 
the  simulation  of  the  device  for  optimum  superlattice  thickness,  the  superlattices  were 
grown  by  molecular  beam  epitaxy  (MBE)  on  semi-insulating  GaAs  substrates.  The 
reactor  is  an  Intevac  Modular  Gen  II  MBE  machine  with  uncracked  As  and  Sb,  and 
elemental  Ga,  In,  and  A1  source  material.  A  4  pm  GaSb  buffer  layer  was  grown  directly 
on  3”  GaAs  substrates.  The  wafer  was  then  broken  into  ~1  cm2  pieces  and  indium- 
mounted  to  molybdenum  blocks.  InAs  is  found  to  have  a  very  narrow  window  for  planar 
growth,  while  high  quality  GaSb  can  be  grown  in  a  wider  range  of  growth  conditions 
when  reflection  high  energy  electron  diffraction  (RHEED)  showed  a  1x3  reconstruction 
pattern  (see  Figure  9).  The  optimum  growth  conditions  for  InAs  layers  were  found  to  be: 
T=400°C  according  to  a  pyrometer,  a  V  to  III  incorporation  rate  ratio  «  3,  and  a  growth 


rate  of  0.5  monolayer/s.  In  this  condition,  RHEED  showed  2x4  reconstruction  patterns. 
The  pyrometer  is  calibrated  with  the  temperature  of  the  transition  from  a  1x5  to  a  1x3 
reconstruction  pattern  in  the  GaSb  buffer  layers.  Based  on  the  theoretical  modeling  and 
simulation  we  chose  the  optimum  structure  for  a  room  temperature  detector  at  X—  1 1  pm. 
The  structure  consisted  of  a  t  =  2  pm  superlattice  with  48  A  InAs,  30  A  GaSb  ,one 
monolayer  of  InSb  at  the  interfaces  as  it  is  shown  to  improve  the  optical  and  electrical 
quality  of  the  superlattice20.  Finally,  the  superlattice  was  capped  with  a  thin  200  A  GaSb 


Figure  9.  The  RHEED  of  the  GaSb  surface  shows  the  lx  pattern  (a),  and  the  3x  pattern 
(b)  these  patterns  are  90°  apart  in  the  azimutal  direction. 

II.G.  Characterization 

The  quality  of  the  superlattices  was  assessed  by  structural,  electrical,  and  optical 
characterization. 

II.G.1.  Structural  Characterization 

High  resolution  x-ray  diffraction  was  used  to  investigate  the  structure  of  the 
material.  X-ray  diffraction  simulation  has  also  been  performed  to  verify  the  superlattice 
structures.  Figure  10  shows  the  typical  x-ray  diffraction  of  the  Type  II  superlattices  and 
the  simulation  results.  Sharp  x-ray  peaks  up  to  the  5th  order  satellite  were  observed. 

This  is  an  indication  of  high  crystal  quality  and  superlattice  uniformity  in  the  samples, 
despite  more  than  6%  mismatch  between  the  lattice  constant  of  the  superlattices  and  the 
GaAs  substrates21.  This  figure  also  shows  the  reproducibility  of  the  growth,  since  the 
second  sample  was  grown  a  week  after  the  first  one. 


Figure  10.  High  resolution  x-ray  diffraction  of  two  samples  and  the  simulation  result. 
Although  one  of  the  samples  was  grown  one  week  after  the  other,  they  are  almost 
identical  and  the  simulation  is  also  in  good  agreement  with  them. 


Figure  10  shows  that  the  mismatch  between  the  average  lattice  constant  of  the  SL 
and  the  lattice  constant  of  the  GaSb  buffer  layer  is  about  0.3%.  The  full  width  at  half 
maximum  (FWHM)  of  the  peeks  is  about  180  arcsec.  Although  high  quality  SLs  with 
thickness  of  about  0.5  pm  can  be  grown  with  such  mismatch,  one  has  to  reduce  the 
mismatch  for  thicker  superlattices.  The  average  lattice  constant  of  the  superlattice  can  be 
adjusted  by  changing  the  thickness  of  the  InAs  layer,  however  this  will  also  change  the 
bandgap  of  the  SL!  Interface  modification  provides  a  better  method  to  change  the 
average  lattice  constant  of  the  superlattices  with  a  little  effect  on  their  bandgap.  As  was 
mention  before,  the  interfaces  can  be  either  InSb-type  or  GaAs-type.  Although  these 
layer  are  very  thin  (about  1  monolayer),  they  have  a  considerable  effect  on  the  average 
lattice  constant  of  the  SL,  since  the  mismatch  between  InSb  and  GaSb  is  about  +0.6% 
and  the  mismatch  between  GaAs  and  GaSb  is  -0.6%.  Considering  the  thickness  of  the 
interfaces,  a  precise  adjustment  of  these  layers  is  obviously  a  difficult  task.  We  managed 
to  routinely  grow  SLs  with  very  low  mismatch  to  the  GaSb  buffer  layers  (see  Figure  1 1) 
by  utilizing  two  techniques.  The  atomic  beam  fluxes  Ga,  In,  Sb,  and  As  were  calibrated 


with  dynamic  RHEED  technique  to  below  1%  error  before  each  growth.  Also,  a  low 
growth  rate  of  about  0.5  monolayer  per  second  was  used,  and  therefore  the  deposition 
time  of  the  interface  was  about  two  seconds  which  is  considerably  longer  than  the 
actuation  time  of  the  shutters  in  our  MBE  system(  ~0.1  second). 


Figure  11.  High  resolution  x-ray  diffraction  of  a  superlattice  grown  on  a  200  A  AlSb 
layer  on  a  GaSb  buffer  layer.  The  SL  has  a  lattice  mismatch  of  less  than  0.05%  to  the 
GaSb  buffer  layer. 

The  surface  morphology  of  the  samples  was  also  assessed  with  a  Digital 
Instruments  Nanoscope  Ilia  atomic  force  microscope.  The  root  mean  square  (rms) 
roughness  of  the  surface  is  less  than  5  A  (see  Figure  12)  which  is  comparable  to  the  rms 
roughness  of  the  samples  grown  on  lattice  matched  GaSb  substrates. 


Figure  12.  (a)  AFM  image  of  the  surface  of  a  GaSb  buffer  on  GaAs  substrate  shows  a 
root  mean  square  (rms)  roughness  of  about  4.2  A.  (b)  The  surface  morphology  of  the 
sample  after  the  growth  of  the  superlattice  on  top  of  the  buffer  layer.  The  rms  roughness 
is  about  4.9  A. 


1I.G.2.  Optical  Characterization 

A  Mattson  Instrument  Galaxy  3000  Fourier  transform  infrared  (FTIR) 
spectrometer  was  used  to  obtain  the  optical  characteristics  of  the  superlattices.  Figure  13 
shows  the  room  temperature  optical  absorption  spectra  of  the  superlattice22.  The  effect  of 
the  substrate  and  GaSb  buffer  layer  was  removed  by  measuring  the  background  with  a 
substrate  and  GaSb  buffer  layer.  The  oscillatory  features  between  0.3  eV  to  0.5  eV  are 
due  to  the  interference  of  internal  reflections  inside  the  SL  (Fabry  Perot  oscillation). 

Since  the  light  path  length  during  the  background  measurement  was  not  identical  to  the 
sample  measurement,  we  could  not  eliminate  the  atmospheric  absorption  effects.  The 
features  around  0.175  eV  (~  6-7  pm)  are  artifacts  due  to  the  high  atmospheric  absorption 
in  this  wavelength  from  water  vapor  absorption  lines. 


Figure  13.  Room  temperature  absorption  spectra  of  a  0.5  pm  thick  sample.  The 
absorption  at  0.12  eV  is  comparable  to  the  absorption  of  HgCdTe  with  similar  bandgap. 


II.G.3.  Electrical  Characterization 


Since  electrons  and  holes  are  confined  in  InAs  and  GaSb  respectively,  we  could 
not  use  the  Hall  measurement  technique  on  superlattices  because  of  the  high  sheet  density 
of  electrons  and  holes  in  the  InAs  and  GaSb  layers.  The  overall  Hall  coefficient  for  this 
material  Rh,  and  conductance  p,  with  two  dominant  channels  of  electrons  and  holes  can 


be  approximated  as: 
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B  is  the  magnetic  field  density  in  the  Hall  measurement,  and  pe  and  (ih  are  the  electron 
and  hole  mobilities  and  n  and  p  are  their  concentrations.  The  formula  shows  clearly  that 
the  overall  Hall  mobility  can  be  much  smaller  than  the  real  mobility  of  electrons  or  holes 
if  <jxy  approaches  zero.  As  the  transport  of  electrons  is  more  important  for  photodetector 
operation  due  to  its  higher  mobility,  single  quantum  wells  of  InAs  were  grown  to 
examine  this  property.  This  provided  a  simple  and  fast  method  for  electrical  assessment 
of  the  interfacial  quality.  After  the  optimization  of  growth  conditions,  the  in-plane  room 
temperature  mobility  of  electrons  in  a  75  A  InAs  well  increased  from  5,000  to  14,000 
cm2/v-s  which  is  about  the  half  the  value  of  bulk  InAs. 

For  superlattices,  field-dependent  Hall  measurement  was  used  ( in  collaboration 
with  Wright  Laboratory,  Wright  Patterson  Air  Force  Base)  to  extract  the  mobility  and 
concentration  of  electrons  and  holes  at  different  temperatures. 

Hall  coefficient  RH,  and  lateral  conductance  p,  were  measured  at  eight  different 
magnetic  flux  densities  for  a  Type  II  superlattice  with  InAs/GaSb  (48A/48A).  A 
minimum  root  mean  square  method  was  used  for  fitting  with  n,  p,  p.e>  Ph>  as  the  fitting 
parameters  to  the  measured  data.  Fitting  is  an  accurate  method,  since  there  are  only  four 
unknowns  for  sixteen  equations.  Figure  14  shows  the  results  of  this  fitting  at  T=300K. 
The  mobility  and  concentration  of  the  electrons  were  calculated  as:  pn~4550  cm  /v.s  and 
n~1.29xl017cm"3  and  for  holes:  pp~1540  cnfVv.s  and  p~1.2x!017  cm'3  at  this  temperature. 


Magnetic  flux  density  (T) 


Figure  14.  The  measured  Hall  coefficient  Rh,  and  conductance  p  of  the  superlattice 
versus  B,  and  the  calculated  values  based  on  the  parameters  in  the  inset  box. 

Figure  15  shows  the  results  of  similar  calculation  at  T  =  120  K.  The  fitted  curve 
clearly  deviates  from  the  measured  p  at  B~1 .25  T.  This  is  due  to  the  fact  that  our  model 
only  considers  electrons  and  holes  with  unique  mobilities,  while  in  reality  each  carrier 
can  have  several  channels  with  different  mobilities. 


Magnetic  flux  density  (T) 


Figure  15.  The  measured  Flail  coefficient  RH,  and  conductance  p  of  the  superlattice 
versus  B  at  T=120  K.  The  fitted  curve  deviates  slightly  from  the  measured  values  of  p  at 
B-1.25  T. 

Figure  16  shows  the  mobility  "spectrum"  of  the  carriers  in  the  superlattice, 
calculated  with  a  program  developed  at  Wright  Laboratory.  It  shows  the  concentration  of 
the  different  carriers  versus  their  mobilities.  Besides  the  peak  at  p^O  which  is  due  to  an 
artifact,  three  other  peaks  are  visible  that  correspond  to  an  electron-type  carrier  with 
mobility  of  about  -5000  cm2/v.s  and  concentration  of  about  ~5xl015  cm'3,  and  two  hole- 
type  carriers  with  mobilities  of -3200  and  -15000,  and  concentrations  of  2xl016cm'3, 
and  4xl014  cm'3  respectively. 


Figure  16.  The  calculated  carrier  concentration  versus  mobility  based  on  the  field- 
dependent  Hall  measurement  method  at  T=120K. 

II.H.  Device  Processing 


II.H.l.  Cleaning 

The  first  step  in  the  processing  is  the  cleaning  of  the  sample  surfaces  from 
mineral-based  and  organic-based  contaminants.  This  is  an  essential  step  of  the 
processing,  since  such  contamination  can  act  as  a  mask  during  chemical  etching  and  also 
prevent  maximum  adhesion  of  the  metal  contacts.  We  noticed  that  a  commonly  used 
degreasing  solution,  Trichloroethylene,  attacks  the  surface  of  samples  that  are  terminated 
with  GaSb.  Therefore,  the  following  cleaning  recipe  was  used: 

1-  Rinse  in  acetone  for  1  min. 

2-  Rinse  in  boiling  methanol  for  1  min. 

3-  Rinse  under  DI  water  for  3  min. 

4-  Repeat  steps  1-3 

5-  Blow-dry  the  samples  with  dry  nitrogen. 


1I.H.2.  Photolithography 

The  first  step  of  the  lithography  process  is  the  photoresist  (PR)  spinning,  in  which 
a  thin  layer  of  PR  is  formed  on  the  surface.  The  thickness  of  the  PR  is  determined  from 
the  PR  viscosity  and  the  rotation  speed  and  time.  We  used  an  AZ5214PR  and  rotation 
speed  of 4000  rpm  to  have  a  PR  thickness  of  about  1 .5pm.  The  next  step  is  a  short 
baking,  known  as  soft  baking,  of  the  samples  at  about  96°C  for  50  seconds.  The  PR 
coated  sample  is  then  exposed  with  UV  light  through  a  chrome-patterned  glass  mask  for  a 
few  seconds.  There  are  two  types  of  PR  processes:  a  positive  process  that  makes  a 
duplicate  of  the  original  mask  and  a  negative  process  that  makes  a  negative  pattern  of  the 
mask.  For  the  positive  process,  the  samples  are  dipped  into  a  developer  solution  and  the 
exposed  areas  will  be  etched  away.  For  a  negative  process  (liftoff),  the  samples  are 
baked  at  about  1 16°C  for  50  seconds,  exposed  to  UV  for  about  1  minute  (without  any 
mask),  and  then  dipped  into  the  developer.  The  last  step  for  a  positive  process  is  a  post 
baking  at  about  1 16°C  for  50  seconds  which  improve  the  adhesion  of  the  PR.  An 
important  difference  between  the  positive  and  negative  process  is  the  shape  of  the  PR  at 
the  edge  of  the  exposed  area.  A  negative  process  makes  an  mushroom  shape  that  is 
necessary  for  a  successful  liftoff  process  which  is  mainly  used  for  the  definition  of  the 
Ohmic  contacts  (see  Figure  17Figure  17.) 


Figure  17.  The  major  processing  steps  for  the  definition  of  the  mesa  (positive  process) 
and  metal  contacts  (negative  process). 


II.H.3.  Etching 

The  grown  material  was  etched  into  mesas  with  wet  (chemical)  and  dry  etching. 
Several  chemical  solutions  that  are  commonly  used  for  etching  of  both  GaSb  and  InAs 
were  studied  for  uniform  and  controllable  etching  of  the  InAs/GaSb  superlattices.  These 
included  Br2:Methanol,  HF:HN03:CH3-C00H,  and  CH3-C00H:HCl:HN03:Br2. 
However,  H3P04:H202:H20  (1:1:10),  which  is  used  for  the  etching  of  InAs,  showed  the 
best  etching  properties.  This  solution  also  provides  a  convenient  etching  rate  of  about 
0.6pm/min.  For  InAs/GalnSb  superlattices  a  better  etching  solution  is  based  on  C6Hg07 
(Citric  acid )  and  H202.  Citric  acid  is  first  dissolved  in  water  with  a  1  gram/ml  ratio,  the 
solution  is  then  mixed  with  an  equal  volume  of  H202  and  stored  for  a  few  minutes  before 
etching.  While  the  etching  rate  of  about  300A/s  is  relatively  slow,  this  solution  provides 
a  very  controllable  etching. 


(a)  (b) 

Figure  18.  (a)  SEM  cross  section  image  of  a  superlattice  etched  with  CeH^:  H2O: 
H202,  and  (b)  High  resolution  SEM  image  of  the  edge  of  the  mesa  shows  a  clean  and 
smooth  etching. 

Dry  etching  techniques  were  also  developed  for  a  higher  uniformity  over  large 
areas,  and  more  vertical  sidewalls.  Several  different  methane  and  chlorine  based  recipes 
were  studied.  A  mixture  of  BCfyAr  (10  sccm:40  seem)  at  a  pressure  of  about  lmTorr, 

RF  power  of  lOOWatt,  and  microwave  power  of  about  800  Watt  showed  the  best 
performance.  It  provided  vertical  sidewalls  as  well  as  (see  Figure  19)  smooth  surfaces 
under  the  etched  areas.  Figure  20  (a)  and  (b)  compares  the  AFM  image  and  surface 
roughness  of  a  superlattice  before  and  after  dry  etching.  The  roughness  is  essentially  the 


Figure  19.  SEM  image  of  the  sidewalls  of  a  mesa  after  ~  lpm  dry  etching. 


(a)  (b) 


Figure  20.  AFM  image  of  the  surface  of  a  superlattice  (a)  before  etching  rms  surface 
roughness  is  7.8  A,  and  (b)  after  ~1  pm  dry  etching  the  rms  surface  roughness  is  9.5  A. 

After  dry  etching  the  photoresist  was  removed  with  a  special  chemical,  AZ  400T 
photoresist  stripper,  provided  by  the  photoresist  supplier.  Although  SEM  and  optical 
microscope  shows  a  clean  surface  after  this  process,  AFM  image  of  the  surface  (see 
Figure  21 -a)  shows  that  there  are  still  a  considerable  photoresist  residuals.  This  tiny  layer 
can  reduce  the  performance  of  the  device,  since  it  prevents  a  good  adhesion  between  the 
metal  contact  and  the  superlattice  in  the  metallization  step.  An  oxygen  ashing  process 
was  used  to  remove  the  residual  photoresist.  It  uses  an  oxygen  flow  rate  of  20  seem  and 
pressure  of  50  mTorr  with  a  RF  power  of  100  Watts  for  40  seconds.  Figure  21-b  shows 
the  AFM  image  of  the  surface  after  the  ashing  process.  The  surface  roughness  is 
comparable  to  the  roughness  of  the  superlattice  before  any  process  (see  Figure  20-a). 
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Figure  21.  AFM  image  of  the  surface  of  the  superlattice  (a)  after  the  photoresist  is 
removed  with  stripper  the  rms  surface  roughness  is  31.5  A,  and  (b)  rms  surface  roughness 
is  7.9  A  after  the  oxygen  ashing  process. 


II.H.4.  Metallization 

A  standard  liftoff  technique  was  used  for  the  definition  of  the  top  metal  contacts 
as  described  in  the  lithography  section.  Ti/Au  (400A/1200A)  were  evaporated  with  an 
electron  beam  evaporator  for  the  contacts.  Figure  22  shows  the  schematic  diagram  and 
the  SEM  image  of  the  processed  device. 


Figure  22.  (a)  Schematic  digram,  and  (b)  SEM  image  of  the  processed  device. 

II.I.  Device  Measurements 

The  samples  were  then  glued  to  a  copper  heatsink  with  an  epoxy  and  attached  to  the 
cold  finger  of  a  liquid  nitrogen  cryostat  equipped  with  a  temperature  controller. 

1I.I.1.  Responsivity 

Spectral  photoresponse  of  the  device  was  measured  using  a  Galaxy  3000  FTIR 
spectrometer  system.  The  samples  were  illuminated  through  the  front  side  with  normal 
incidence.  The  absolute  response  of  the  photodetectors  was  calculated  using  a  blackbody 
test  set,  which  is  composed  of  a  blackbody  source  (Mikron  305),  preamplifier  (EG&G 
PA-6),  lock-in  amplifier  (EG&G  5209),  and  chopper  system  (Stanford  Research  System 
SR540).  Figure  23  shows  the  responsivity  of  the  device  in  the  2-  17pm  wavelength  range 
at  78K  and  300K  with  an  in-plane  electrical  field  of  5V/cm. 


Figure  23.  The  responsivity  spectra  of  the  device  at  78K  and  300K  with  an  in-plane 
electrical  field  of  5V/cm. 

To  assess  the  temperature  dependence,  the  current  responsivity  of  the  device  was 
measured  at  AM  0.6pm  wavelength  from  78K  to  room  temperature  at  a  constant  electrical 
field.  Figure  24  shows  the  responsivity  of  the  detector  at  A,=l  0.6  pm  versus  the  detector 
temperature.  In  order  to  see  whether  current  responsivity  follows  a  power  function,  we 
fit  the  data  to  an  allometric  function.  An  allometric  fit  has  a  general  form  of  y=AxB  where 
x  is  the  variable  and  A  and  B  are  the  fitting  parameters.  This  fit  shows  that  the 
responsivity  of  the  detector  is  nearly  proportional  to  T'  . 
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Figure  24.  The  current  responsivity  of  the  device  versus  temperature  at  X,=l  0.6pm  at 
constant  voltage  bias.  The  squares  are  the  measured  points  and  the  line  is  an  Allometric 
fit  (AT8)  to  the  points  which  shows  that  responsivity  is  nearly  proportional  to  T  . 

This  is  an  unusual  behavior  since  responsivity  of  the  narrow  gap  material  is  usually 
an  exponential  function  of  temperature  at  higher  temperatures  where  Auger  ^ 

recombination  is  the  dominant  recombination  mechanism.  Theoretically,  Zegrya  et  al 
showed  that  the  Auger  recombination  rate  is  a  power  function  of  the  temperature 
(proportional  to  T2)  in  Type  II  heterostructures  compared  to  the  exponential  function  in 
the  bulk  semiconductors.  Since  the  current  responsivity  is  proportional  to  the  carrier 
lifetime  (which  is  dictated  by  Auger  recombination  rate,  this  power  dependency  of  the 

responsivity  indicates  a  good  agreement  with  the  theoretical  predictions. 

The  effective  lifetime  of  the  carriers  was  also  extracted  from  the  responsivity  and 
Hall  measurements  on  a  t=0.5pm  thick  superlattice: 

t  = _ !i _  (27) 

ef  E(pe+pp) 

where  l=2mm  is  the  device  length,  g  is  the  photoconductor  gain,  and  E=5V/cm  is  the 
electrical  field.  The  gain  of  the  device  can  be  calculated  from: 


(28) 


R:hc 

g  =  -T - 

where  Rj=2mA/W  is  the  current  responsivity,  h  is  the  Planck  constant,  c  is  the  speed  of 
light,  X=10.6pm  is  the  wavelength  of  the  light,  rj  is  the  quantum  efficiency,  and  q  is  the 
electron  charge.  Assuming  an  internal  quantum  efficiency  near  unity  and  negligible 
reflection  from  the  bottom  of  the  superlattice  and  unpolished  backside  of  the  substrate, 
the  quantum  efficiency  can  be  calculated  from: 

r)  =  (1  —  r)(l  -  e”“‘)  (29) 

where  r  is  the  top  surface  reflection  coefficient,  a=1.8xl03cm  1  is  the  absorption 
coefficient  of  the  superlattice  and  t-  0.5  um  is  the  thickness  of  the  superlattice.  Assuming 
r»0.3,  the  quantum  efficiency,  photoconductive  gain,  and  carrier  lifetime  can  be 
calculated  from  above  formulas  as:  r|=6.02%,  g=3.9xl0 3,  and  Tef=26.8nsec.  The  effective 
lifetime  is  about  an  order  of  magnitude  longer  than  the  carrier  lifetime  in  HgCdTe  ^ 
photoconductors  with  similar  bandgap  and  carrier  concentration  at  room  temperature  . 
Since  Auger  recombination  is  the  dominant  recombination  mechanism  at  room 
temperature,  we  believe  that  the  enhancement  of  carrier  lifetime  is  due  to  the  suppression 
of  Auger  recombination  in  the  Type  II  superlattice. 


11.1.2.  Noise  Measurements 

The  noise  was  measured  with  a  fast  fourier  transform  (FFT)  spectrum  analyzer 
(Stanford  Research  System  SR  760)  and  a  low  noise,  wide  band  pre-amplifier  with  54dB 
voltage  gain  (EG&G  PA- 100).  Figure  25  shows  the  input  noise  spectrum  of  the  FFT 
analyzer,  the  output  noise  spectrum  of  the  shorted  pre-amplifier,  and  the  output  noise  of 
the  pre-amplifier,  when  it  was  connected  to  the  detector.  The  detector  was  biased  by  the 
pre-amplifier  at  Vb=5  volts. 


Figure  25.  The  measured  input  noise  spectrum  of  the  fast  Fourier  transform  (FFT) 
spectrum  analyzer,  the  output  noise  spectra  of  the  pre-amplifier  with  a  shorted  input  and 
the  output  noise  spectra  of  the  pre-amplifier  when  it  is  connected  to  the  detector  at  300K. 


The  mean-square  noise  of  the  detector  can  be  modeled  as: 


V2  =  V2  +  V?  +  V, 


GR  J 


(30) 


where  Vn  is  the  overall  noise  of  the  detector,  Vj  is  the  Johnson-Nyquist  noise,  Vi/f  is  the 
1/f  noise  and  Vqr  is  the  generation-recombination  noise.  The  value  of  the  Johnson  noise 
can  be  calculated  as 

V2  =  4kTRAf ,  (31) 

where  k  is  the  Boltzmann  constant,  T  is  the  temperature  and  R=76  ohms  is  the  resistance 
of  the  device.  The  value  of  Johnson  noise  for  the  device  at  room  temperature  is 
1.12nV/Hz1/2.  The  1/f  noise  can  be  approximated  with: 


(32) 


where  fi/f  is  a  constant  which  depends  on  the  sample,  and  f  is  the  frequency.  This  shows 
that  at  high  enough  frequencies,  1/f  noise  can  be  negligible  compared  to  the  other  two 


types  of  noise.  Then  the  value  of  the  generation-recombination  noise  can  be  extracted 
from  the  total  noise  of  the  device  in  this  range  and  the  value  of  the  Johnson  noise.  The 
noise  equivalent  circuit  model25  was  used  to  extract  the  noise  of  the  photodetector  as 
Vn=1.7  nV/Hz1/2  above  10kHz.  From  the  above  equations,  the  value  of  the  generation- 
recombination  noise  can  be  calculated  as  Vgr=1  .28nV/Hzl/2.  The  generation- 
recombination  noise  can  be  approximated  as: 

v  2Vb  1  +  b  np  xAf  ]  (33) 

GR  (lwt)/2  bn  +  p  l^n  +  p  1  +  co2t2J 

where  w=4mm  is  the  detector  width,  co  is  the  angular  frequency  and  ot«1  in  this 
experiment,  and  b  is  the  ratio  of  |_ic  to  Up  .  From  the  Vgr  equation,  the  value  of  the  carrier 
lifetime  can  be  calculated  as  x=17nsec  which  is  close  to  the  value  of  the  lifetime 
extracted  from  the  optical  response  measurement.  It  should  be  noted  that  the  calculation 
of  the  carrier  lifetime  from  the  electrical  noise  of  the  device  is  not  an  accurate  technique 
and  may  only  provide  a  rough  estimation  for  the  carrier  lifetime. 


ILL 3.  Response  Time  Measurement 


As  was  explained  in  section  V.A.,  one  of  the  motivations  for  the  realization  of 
uncooled  photon  detectors  is  the  need  for  uncooled  detectors  with  response  times 
x  <  30  msec.  Therefore,  it  was  very  important  to  measure  the  response  time  of  the 
Type  II  uncooled  detectors.  However,  conventional  methods  such  as  using  mechanical 
choppers  could  not  achieve  time  accuracy  below  the  millisecond  range.  We  used  a  room 
temperature  quantum  cascade  laser,  developed  at  CQD26,  as  a  narrowband  and  high  speed 
infrared  source  in  our  measurement.  The  schematic  diagram  of  the  setup  is  shown  in 
Figure  26.  The  pulse  generator  and  laser  driver  were  inside  an  Avtech  AVR-4A-PW 
which  is  capable  of  generating  high  power  electrical  pulses  with  fall  time  of  about  5  nsec. 
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Figure  26.  The  schematic  diagram  of  the  time  response  measurement  setup. 


The  quantum  cascade  laser  was  uncooled  and  operated  at  X— 8. 5 pm  with  a  negligible  time 
delay.  An  EG&G  PA- 100  low-noise  pre-amp  was  used  to  amplify  the  detector  signal. 
Unfortunately,  the  pre-amp  is  not  very  fast  and  has  a  fall  time  of  more  than  40  nsec.  The 
output  signal  of  the  pre-amp  was  measured  with  a  Tektronix  TDS  520B  digital 
oscilloscope.  It  shows  a  90%  to  10%  fall  time  of  about  68  nsec,  as  shown  in  Figure  27, 
for  the  whole  setup  and  hence  the  detector  has  a  fall  time  below  40nsec  . 


Figure  27.  The  overall  response  of  the  setup. 


II.  J.  Device  Performance  and  Comparison  with  State  of  the  Art 

The  performance  of  infrared  detectors  is  usually  compared  based  on  their 
detectivity  which  is  an  indication  of  their  sensitivity.  Knowing  the  responsivity  and  the 
noise  of  the  devices,  we  could  calculate  their  detectivity.  The  uncooled  devices  show  a 
detectivity  of  about  1 ,3xl08cmHz1/2/W  at  X=llpm  which  is  higher  than  the  detectivity  of 
commercially  available  HgCdTe  at  similar  wavelengths  and  temperature27  (see  Figure 
28).  Unlike  HgCdTe,  these  Type  II  superlattices  are  grown  on  conventional  GaAs 
substrates,  and  hence  highly  uniform  material  can  be  grown  on  three  and  five  inch  wafers 
readily. 


Calculated  Auger  Limited  MCT  Detectivity  with 
MCT  and  Type-II  superlattice  experimental  result 
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Figure  28.  The  detectiviy  of  Type  II  superlattices  compared  to  the  theoretical  limit  and 
experimental  detectivity  of  HgCdTe  (MCT)  detectors  at  300  K  at  A,~l  1  pm. 


In  comparison  to  the  thermal  detectors,  such  as  microbolometers,  Type  II 
superlattices  have  similar  detectivity,  but  are  showing  at  least  five  orders  of  magnitude 
faster  response28.  Although  thermal  detectors  with  higher  detectivity  are  possible,  the 
price  that  one  has  to  pay  is  the  speed.  The  change  of  the  temperature  of  a  bolometer  is: 

AT  =  (34) 

Gyl  +  co2x~ 

where  AT  is  the  temperature  change,  P0  is  the  IR  power  that  comes  to  the  bolometer 
surface,  r|  is  the  percentage  of  the  IR  power  that  is  absorbed,  G  is  the  thermal 
conductance  of  the  bolometer,  co  is  the  frequency  of  the  IR  emission,  and  x  is  the  thermal 
response  time  of  the  bolometer  and  is  x=C/G,  where  C  is  the  thermal  capacitance  of  the 
bolometer.  It  is  easy  to  find  that  at  low  frequencies  cox«l ,  the  change  of  the  temperature 
is  AT~r)P0/G.  Therefore,  the  major  effort  to  improve  the  responsivity  of  the  thermal 
detectors  is  based  on  the  reduction  of  the  thermal  conductance,  G.  However,  this  will 
increase  the  response  time  of  the  device  if  one  can  not  reduce  its  thermal  capacitance, 
since  x=C/G.  Unfortunately,  the  thermal  capacitance  of  the  microbolometers  cannot  be 
reduced  further  since  it  requires  a  thinner  or  smaller  device  that  is  not  practical.  In  order 
to  decrease  the  thermal  conductance  of  the  device,  the  length  of  the  legs  that  are 
supporting  the  sensitive  layer  can  be  increased  (see  Figure  29). 
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Figure  29.  Schematic  diagram  of  microbolometers  with  (a)  orginal  pixel  legs,  and  (b) 
two  times  longer  legs. 

Figure  30  shows  the  relative  frequency  response  of  pixels  with  different  leg  lengths 
presented  by  Raytheon.  Although  increasing  the  leg  length  by  a  factor  of  four  decreases 
the  NEDT  by  nearly  a  factor  of  two,  the  3  dB  frequency  knee  also  decreases  by  a  factor 
of  four29. 
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Figure  30.  The  relative  frequency  response  of  different  test  pixels  with  standard 
Raytheon  Sb-151  leg  length  and  2  times  and  4  times  longer  legs. 


III.  Realization  of  Uncooled  Type  II  Photodiodes  in  the  Long 
Wavelength  IR  (LWIR)  range 

III.A.  Motivation 

The  need  for  uncooled  photon  detectors  in  the  mid  and  long  wavelength  IR  ranges 
as  well  as  the  limitation  of  currently  available  photon  detectors  have  been  detailed  in 
section  V.A.  Although  we  demonstrated  a  high  performance  uncooled  photoconductor  in 
the  long  wavelength  range  based  on  Type  II  superlattices,  this  device  is  more  suitable  for 
single-element  detectors  and  cannot  be  easily  used  for  staring  two-dimensional  (2D)  focal 
plane  arrays  (FPA).  Unlike  a  photoconductor,  the  current  of  a  photodiode  flows 
perpendicular  to  the  surface,  and  hence  it  is  scaled  by  the  area  of  the  device.  For 
example,  the  current  of  a  50  pm  x50  pm  photoconductor  at  a  given  voltage  bias  is  equal 
to  the  current  of  a  1  cm  x  1  cm  device,  while  the  current  of  a  50  pm  x50  pm  photodiode 
is  40000  times  less  than  the  current  of  a  1  cm  x  1  cm  device.  This  fact  becomes  very 
crucial  for  currently  used  large  area  FPAs,  with,  for  example,  512  x  512=  262,144 
elements. 

Moreover,  photodiodes  can  operate  even  at  zero  bias  which  not  only  reduces  the 
bias  and  heat  dissipation  requirements  significantly,  but  also  eliminates  the  1/f  noise. 

This  type  of  noise  increases  inversely  with  frequency,  and  can  be  the  dominant  source  of 
the  total  noise  for  low  frequency  applications  such  as  IR  imaging  systems. 

III.B.  Modeling  and  Simulation  of  the  Device 

The  active  layer  of  these  p-i-n  photodiodes  is  a  short  period  superlattice  and  hence 
the  precise  modeling  of  the  device  requires  the  modeling  of  the  carrier  transport  in  such 
superlattices.  However,  such  modeling  is  very  complex  and  requires  massive  numerical 
calculations30.  Therefore,  for  these  superlattice-based  devices,  we  decided  to  adopt  the 
general  models  available  for  bulk  semiconductors.  Fortunately,  the  simulated  results 
based  on  such  modeling  showed  an  excellent  agreement  with  our  experimental 
measurements.  This  means  that  the  models  are  suitable  for  devices  with  a  wide  range  of 
cutoff  wavelengths  (from  ~8  pm  to  ~  1 6  pm)  as  well  as  operating  temperatures  (80  K  to 
300  K).  The  total  current  density  of  a  photodiode,  Jtotai,  is  the  summation  of  the  dark 
current  density,  Jdark,  which  is  a  function  of  the  voltage  bias  and  the  photocurrent  density, 
Jphoto,  which  is  a  function  of  the  bias  and  photon  flux.  We  simulated  photo  and  dark 
current  separately,  since  they  can  be  measured  independently. 


III.B.  1.  Dark  Current  Modeling(Voltage-Current  Characteristics) 

The  main  current  components  of  a  photodiode  under  no  illumination  (dark 
current)  are  the  diffusion,  the  generation-recombination  (GR),  and  the  tunneling  currents: 

Lark  =Liff  Lr  L  @5) 

In  the  following,  we  develop  the  formalism  needed  for  the  modeling  of  the  current 
density  of  these  components. 


a,  Diffusion  Cu[rent  Density 

The  diffusion  current  is  a  fundamental  mechanism  in  a  p-n  junction  and  is  due  to 
the  diffusion  of  the  minority  carriers  to  edge  of  the  depleted  layer.  The  density  of  the 
diffusion  current  Joiff  is  related  to  the  voltage  bias  of  the  diode  by: 


where  VB  is  the  voltage  bias,  Js  is  the  leakage  current  density  of  the  diode,  q  is  the 
electron  charge,  k  is  the  Boltzmann  constant,  and  T  is  the  absolute  temperature.  The 
leakage  current  density  can  be  calculated  from31: 


qD  nf  , 
+  -  tanh 


qDhnf 


where  n  and  p  are  the  electron  and  hole  concentrations  in  the  p-type  and  n-type  layers 
respectively,  De  and  Dh  are  their  diffusion  coefficients,  and  Le  and  Lb  are  their  diffusion 
lengths.  xn  is  the  thickness  of  the  n-type  layer  and  xp  is  the  thickness  of  the  p-type  layer, 
nj  is  the  intrinsic  carrier  concentration  and  is  calculated  based  on  a  three  dimensional 
density  of  states  approximation: 
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where  me  and  mh  are  the  electron  and  hole  effective  masses.  Eg  is  the  effective  bandgap 
of  the  superlattices.  The  diffusion  lengths  can  be  extracted  from  the  carrier  mobility  p  at 
a  given  temperature  T,  using: 

D-— n  (39) 

q 

Also,  the  diffusion  length  can  be  calculated  from  the  carrier  lifetime  x  and  the  diffusion 
coefficient  D: 

L  =  Vm  (40) 


Generation-Recombination  Current  Density 


The  generation-recombination  (GR)  current  is  due  to  the  GR  process  in  the 
depleted  layer  of  the  photodiode.  Our  approximation  was  based  on  a  symmetric  p-n 
junction32,  since  it  is  much  simpler  than  the  approximation  for  a  n-p+  junction  while  it 


gives  similar  results  for  our  parameters: 


where  tgr  is  the  GR  carrier  lifetime  in  the  depleted  layer.  w(Vb)  is  the  thickness  of  the 
depleted  layer  at  a  voltage  of  VB,  and  using  space  charge  approximation  it  is: 
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where  E(VB)  is  the  electric  field  in  the  space  charge  region  and  can  be  approximated  as: 
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where  8s  is  the  permittivity  of  the  semiconductor. 
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c.  Tunneling.  Current  Density. 

Tunneling  current  is  due  to  the  tunneling  of  the  electrons  from  the  valence  band  to 
the  conduction  band  through  the  bandgap  of  the  semiconductor.  Here,  the  tunneling 
current  of  the  detectors  is  approximated  based  on  a  triangular  potential  barrier  : 
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III.B.2.  Photocurrent  Modeling  ( Quantum  Efficiency) 

The  density  of  the  photocurrent  of  a  photodiode  under  a  photon  flux  density  of  <f> 
is  directly  related  to  its  quantum  efficiency,  q  : 

J Photo  (44) 

where  q  is  the  electron  charge.  The  quantum  efficiency  of  a  p-i-n  photodiode  is  the 
summation  of  the  efficiency  of  the  n-type  layer,  depleted  layer,  and  the  p-type  layer: 

fl=nn  +riw  +'nP  (45) 
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Figure  31.  Schematic  of  a  p-i-n  photodiode,  w  is  the  thickness  of  the  depleted  layer, 
while  xn  and  xp  are  the  thickness  of  the  undepleted  n-type  and  p-type  layers. 


For  our  devices  that  are  illuminated  from  the  n-side  (see  Figure  31),  the  quantum 
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r|w  =  (1  -  r){e-ax"  -e_a(Xn+w)}  (48) 

where  r  is  the  surface  reflectivity,  ot  is  the  optical  absorption  coefficient,  and  Le  and  Fh 
are  the  electron  and  hole  diffusion  lengths. 

III.C.  Device  Design 


There  are  two  important  goals  in  the  design  of  an  uncooled  photodiode  in  the  LWIR 
range.  First,  the  differential-resistance  area  product  (RoA)  must  be  high  enough  to  reduce 
the  thermal  current  noise  and  the  dark  current  of  the  device.  In  fact,  this  has  been  one  of 
the  major  obstacles  for  the  realization  of  practical  IR  imaging  systems  based  on 
photodiodes.  Second,  the  quantum  efficiency  of  the  device  must  be  high  enough  to  have 
a  considerable  signal  level. 

The  value  of  R0A  can  be  easily  calculated  from  the  devices’  dark  current  density: 


R0A 


1 


vB=0 


(49) 


It  is  also  helpful  to  separate  the  contribution  of  different  components  of  the  dark  current 
to  the  RoA.  These  are  namely  the  diffusion  related  (RoA)diff,  the  GR  related  (RoA)gr,  and 


the  tunneling  related  (RoA)t: 
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where  Js ,  tGr,  VB;  w(Vb),  nj,  me,  Eg,  and  ss  are  defined  in  the  previous  section.  Neff  is  the 
effective  doping  density  which  is  equal  to  np/(n+p),  and  h  is  the  Planck  constant. 

Since  these  resistances  are  effectively  in  parallel,  the  overall  RoA  product  is: 

1  1^1  _ 1  (51) 

R0A  (Rtl  A)diff  (R0A)gr  (R0A)t 

This  means  that  RoA  is  similar  to  the  component  with  the  lowest  value.  In  the  case  of 
devices  with  a  cutoff  wavelength  close  to  Xc=8pm  at  room  temperature,  (RoA)oiff  is  much 
lower  than  the  other  components,  and  hence  is  the  most  important  parameter  of  the 
design.  This  parameter  can  only  be  increased  by  decreasing  the  leakage  current  Js.  A 
quick  study  of  the  formula  for  Js  shows  that  one  can  decrease  Js  by:  first  decreasing 
mobility  p,  second  increasing  carrier  lifetime  x,  and  third  decreasing  xn  and  xp  . 

The  mobility  of  the  electrons  and  holes  are  naturally  reduced  by  interface 
scattering  in  Type  II  superlattices  and  the  estimated  mobility  is  on  the  order  of  ~100 
cm2/ Vs  for  electrons  and  ten  times  lower  for  holes35.  This  is  a  definite  advantage  of 
Type  II  SLs  since  these  are  orders  of  magnitude  smaller  than  similar  values  for  InAsSb  or 
HgCdTe  where  electron  mobility  is  around  pe  -100,000  cm2/V s! 

In  the  previous  chapter,  we  showed  that  the  carrier  lifetime  x  of  a  properly 
designed  Type  II  SL  is  nearly  one  order  of  magnitude  longer  than  the  carrier  lifetime  in  a 
narrow  gap  bulk  material  (e.g.  InAsSb  and  HgCdTe). 

The  third  point  indicates  that  xn  and  xp  must  be  reduced.  However,  a  glance  at  the 
formula  of  quantum  efficiency  shows  that  such  action  can  also  reduce  the  quantum 
efficiency  of  the  device.  The  question,  however,  is  which  value  of  xn  and  xp  is  optimum 
in  this  tradeoff.  Using  the  parameters  that  were  extracted  from  the  simulation  of  the  dark 
current  of  the  processed  devices,  I  could  simulate  the  quantum  efficiency  versus  xn  and  xp 
as  is  shown  in  Figure  32.  The  parameters  that  were  used  are  pe=100  cm  /Vs, 
ph=10  cm2/Vs,  oc=2000  cm'1,  w=45  nm,  and  x=20  ns.  Figure  33  shows  the  calculated 
leakage  current  density  Js  for  the  same  set  of  parameters. 


Figure  32.  The  calculated  quantum  efficiency  of  a  device  with  A.c-8pm  at  room 
temperature  versus  the  thickness  of  the  n-type  (xn)  and  p-type  (xp)  layers. 


The  value  of  leakage  current  density  is  clearly  a  more  sensitive  function  of  layer 
thickness.  For  example,  Js  is  nearly  five  times  higher  for  xn=Xp:=0.2  pm  compared  to 
xn=xp=l  .0  pm,  while  t|  is  only  ~3.6  times  higher. 


Xn(nm) 

Figure  33.  The  calculated  leakage  current  density  (in  A/cm2)  of  a  device  with  A,c=8pm  at 
room  temperature  versus  the  thickness  of  the  n-type  (xn)  and  p-type  (xp)  layers. 

The  detectivity  is  proportional  to  quantum  efficiency  and  the  inverse  square  root  of  the 
leakage  current  density  or  equally  to  the  square  root  of  RoA  product  (Jsocl/RoA). 
Therefore,  the  detectivity  of  a  device  with  xn=xp=1.0  pm  is  3.5/  -\[5  =~  1.56  times  higher 
than  a  device  with  xn=xp=0.2  pm  while  its  quantum  efficiency  is  3.5  times  higher. 

So  far  the  thicker  device  in  this  example  has  clear  advantages:  higher  detectivity 
(56%)  and  higher  quantum  efficiency  (3.5  times),  while  its  only  disadvantage  is  lower 
RoA  (5  times).  However,  I  realized  that  the  thickness  is  a  very  important  issue  which  can 
be  used  to  improve  the  performance  of  the  device.  I  developed  a  novel  design  in  which 
several  p-i-n  devices  can  be  grown  on  top  of  each  other  and  act  as  a  series  of  detectors.  A 
detailed  study  of  this  new  device  design,  called  a  multiple  junction  detector  (MJD)  in  this 
thesis,  showed  excellent  properties  compared  to  the  conventional  p-i-n  design  with  a 
similar  thickness.  For  example,  an  MJD  consisting  of  five  p-i-n  junctions  with 
Xn=Xp=0.2  pm  have  an  overall  thickness  and  growth  time  of  the  device  with 
Xn=Xp=1.0  pm.  However,  the  MJD  has  5x5=25  times  higher  RoA  and  V25  /3.5  =  1.43 
times  higher  detectivity,  but  3.5  times  lower  quantum  efficiency  compared  to  a  simple 
p-i-n  with  similar  thickness. 

Fortunately,  high  quantum  efficiency  is  not  as  important  as  high  detectivity  for 
the  imaging  of  room  temperature  objects  (e.g.  medical  thermal  imaging),  since  the  photon 


flux  is  quite  high.  Good  examples  of  high  performance  IR  imaging  arrays  with  only  a 
few  percent  quantum  efficiency  are  cooled  Schottky  barrier  silicide  detectors  and 
uncooled  microbolometers. 

Moreover,  the  above  values  of  quantum  efficiencies  are  the  so-called  internal 
quantum  efficiencies.  The  “external”  quantum  efficiency  is  the  overall  quantum 
efficiency  of  a  system  consisting  of  the  detector  and  a  read-out  circuit.  The  ratio  of  the 
external  to  the  internal  quantum  efficiency  is  equal  to  the  ratio  of  the  current  that  reaches 
the  read-out  circuit  to  the  photo-generated  current.  This  ratio  is  related  to  the  impedance 
of  the  detector  Rd  and  the  read-out  circuit  Rr  as: 

T external  _ _ j; _  (52) 

internal  1  +  Rr/Rd 

This  formula  shows  that  external  quantum  efficiency  approaches  the  internal  quantum 
efficiency  only  if  the  resistance  of  the  detector  is  much  higher  than  the  read-out  circuit. 

In  fact,  the  differential  resistance  of  the  bulk-based  uncooled  photodiodes  is  so  low  that 
the  above  situation  cannot  be  met  and  external  quantum  efficiency  is  much  lower  than  the 
internal  values.  Therefore,  the  MJD  design  of  the  above  example  with  25  times  higher 
differential  resistance  can  have  even  a  higher  overall  quantum  efficiency  than  the  simple 
p-i-n  design,  depending  of  the  impedance  of  the  read-out  circuit. 

III.D.  Growth 

The  structures  were  grown  using  cracked  As  and  Sb  sources.  The  cracking  zone 
temperature  for  both  cells  was  900°C.  The  deposition  rates  of  the  material  were 
calibrated  with  dynamic  RHEED  oscillation  to  within  1%.  First,  1pm  of  GaSb  contact 
layer  doped  with  Be  (NA=lxl018cm"3)  was  grown  on  a  GaSb-p  substrate.  Then  a  stack  of 
five  devices  was  grown,  each  of  which  consisted  of  20  periods  of  p-type  InAs/GaSb:Be 
(39A/40A),  20  periods  of  nominally  undoped  InAs/GaSb  (39A/40A),  and  20  periods  of 
n-type  InAs:Si/GaSb  (39A/40A)  superlattices.  GaSb  layers  in  this  superlattice  had  a 
graded  doping  from  1018cm'3  to  2xl017cm‘3 .  In  the  n-type  superlattice,  InAs  layers  were 
doped  with  Si  (ND=2xl018cm'3).  The  shutter  sequences  were  designed  such  that  both 
interfaces  were  InSb  type.  Finally,  the  growth  was  capped  with  a  0.01  pm  InAs. Si  layer 
(ND=Txl018cnf3).  The  growth  temperature  was  about  520°C  for  the  GaSb  and  395°C  for 
the  superlattices  according  to  a  pyrometer  calibrated  with  the  surface  reconstruction 
transition  temperatures  of  GaSb  (~390°C)  and  InSb  (~380°C). 

III.E.  Structural  Characterization 

High  resolution  x-ray  diffraction  and  atomic  force  microscopy  (AFM)  showed  excellent 
crystal  quality  and  surface  smoothness  of  the  samples.  We  could  routinely  grow  samples 
with  less  than  0.1%  mismatch  to  the  GaSb  substrates  and  less  than  4  A  root  mean  square 
(rms)  surface  roughness.  Figure  34  shows  the  rocking  cure  diffraction  pattern  of  the 
grown  material  in  the  (004)  crystal  direction.  High  order  satellites  with  sharp  peaks 
indicate  a  uniform  crystal  period  and  excellent  crystal  quality  and  sharp  interfaces. 


Figure  34.  High  resolution  x-ray  diffraction  pattern  of  the  grown  material  shows  the 
small  mismatch  of  the  superlattice  to  the  GaSb  substrate  as  well  as  high  order  satellites 
with  sharp  peaks. 


Figure  35  shows  the  atomic  force  microscopy  (AFM)  image  of  the  surface  as  well 
as  the  results  of  a  roughness  analysis  program.  The  rms  surface  roughness  over  a 
20  pm  x  20  pm  area  is  only  1.9A  which  is  the  best  reported  value  for  Type  II 
InAs/Gai  ,xInxSb/InAs  superlattices  grown  on  GaSb  substrates. 
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Figure  35.  AFM  image  of  the  surface  morphology  of  the  sample  as  well  as  the  roughness 
analysis.  The  rms  surface  roughness  over  a  20  pm  x  20  pm  area  is  only  1.9A. 

High  resolution  AFM  images  of  the  surface  of  the  grown  material  show  wide  atomic 
steps  (see  Figure  36)  which  is  an  indication  of  excellent  surface  smoothness  of  the 
samples. 


Figure  36.  High  resolution  AFM  image  of  the  surface  of  the  device  shows  atomic  steps 
as  wide  as  2  pm  indicating  excellent  surface  smoothness  of  the  grown  material. 

III.F.  Device  Processing 

The  samples  were  then  processed  into  400  pm  x  400  pm  mesas  using  standard 
photolithography  and  a  new  wet  etching  solution.  The  major  challenge  for  the 
development  of  the  etching  solution  comes  from  the  fact  that  it  should  etch  a  superlattice 
with  three  different  materials:  InAs,  InSb,  and  GaSb.  Although  one  can  use  a  solution 
which  only  etches  one  material  and  the  others  will  be  pilled  off  due  to  the  under  etching 
effect,  I  realized  this  will  not  provide  a  controllable  etching  method.  The  only  existing 
etching  solutions  that  etch  nearly  all  EH-V  materials  with  a  similar  rate  are  bromine  based 
solutions.  Unfortunately  the  etching  mechanism  is  diffusion-limited,  and  hence  the 
inevitable  non-uniformity  of  the  solution  flow  makes  the  etching  rate  and  the  shape  of  the 
sidewalls  non-uniform.  They  also  attack  the  photoresist  mask,  and  therefore  a  dielectric 
mask  (e.g.  Si02)  must  be  used  which  increases  the  processing  complexity  considerably. 
The  solution  that  I  developed  is  a  mixture  of  citric  acid,  hydrogen  peroxide,  phosphoric 
acid  and  water  (2:3:2:20).  It  etches  all  superlattice  layers  with  a  similar  etching  rate  of 
~0.5  pm/min  which  is  well  in  the  controllable  range.  It  has  also  a  reaction-limited 
mechanism,  and  so  the  etching  rate  is  independent  of  the  agitation  level  of  the  solution 
and  highly  uniform  mesas  can  be  etched  over  a  large  area  (see  Figure  37). 

For  top  and  bottom  contacts  we  used  Ti/Au  (500  A  /  1200  A),  defined  with  electron 
beam  evaporation  and  lift-off  techniques.  No  passivation  or  anti-reflection  coating  was 
used  on  the  surfaces  of  the  samples. 


(a)  (b) 


Figure  37.  The  SEM  images  of  (a)  a  400  pm  x  400  pm  mesa  and  (b)  the  magnifide  edge 
of  the  mesa.  The  wet  etching  provides  smooth  beveled  edges  with  a  ~60°  angle. 


III.G.  Device  Measurement 


III.G.l.  Current-Voltage  Characteristics 

The  current-voltage  (I-V)  characterization  of  the  devices  was  measured  with  a  HP 
41 55 A  parameter  analyzer.  Figure  38  shows  the  current  density  and  differential  resistance 
area  product  (RA)  versus  the  bias  of  a  device.  The  value  of  the  differential  resistance  at 
zero  bias  is  RoA=1.36xlO'2  Qcm2  which  is  more  than  two  orders  of  magnitude  higher 
than  R0A=~10'4  Qcm2  of  HgCdTe  detectors  with  a  similar  bandgap  at  room 
temperature  . 


Figure  38.  The  current  density  and  differential  resistance  area  product  (RA)  of  a  device 
versus  the  voltage  bias. 

Figure  39  compares  the  calculated  current  density  of  a  device  and  the  measured 
values.  Current  density  simulation  uses  the  formalism  given  in  VI.B.l  and  shows 
excellent  agreement  to  the  experimental  data.  The  fit  parameters  are  only  the  carrier 
lifetime  Te=Th=20  nsec,  and  tgr=0.6  nsec.  We  used  pe=  1 00cm2/V s  and  Ph=10  cm2/V s 
which  were  extracted  from  the  experimental  results  in  Gai-xInxSb/InAs  superlattices  with 
a  similar. 
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Figure  39.  The  calculated  current  density  versus  the  voltage  bias  compared  to  the 
measured  values  for  sample  SSMBE  1224  at  T=300K  shows  excellent  agreement  in 
reverse  and  forward  biases. 


III.  G.  2.  Optical  Responsivity 

Absolute  spectral  responsivity  was  calculated  from  the  measured  spectral 
response  of  the  devices,  using  a  Fourier  transform  infrared  (FTIR)  spectroscopy  system, 
and  the  device’s  photoresponse  to  a  calibrated  blackbody  setup.  The  peak  responsivity 
was  Rj=0. 1 4  A/W  at  X=1  (am  leading  to  a  Johnson  noise  limited  detectivity  of 
D*=  1 ,2xl08  cmHz1/2/W  at  room  temperature.  Figure  40  shows  the  detectivity  of  a  device 
versus  IR  wavelength  and  energy.  Although  the  IR  path  length  in  the  atmosphere  was 
only  about  15  cm,  CO2  and  water  vapor  absorption  features  are  visible  in  the  spectrum. 
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Figure  40.  The  Detectivity  of  a  device  versus  the  wavelength  and  energy  of  the  IR 
radiation.  Features  at  4.3  (am  and  6-7  pm  are  due  to  the  C02  and  water  vapor  absorption 
during  the  measurement. 


III.  G.3.  Noise  Measurement 

Although  non-equilibrium  HgCdTe  and  InAsSb  detectors  with  high  differential 
resistance  under  reverse  bias  have  been  demonstrated  recently37,  a  high  1/f  noise 
degrades  the  performance  of  these  devices  below  several  megahertz  by  two  to  three 
orders  of  magnitude.  Consequently,  these  non-equilibrium  devices  cannot  be  used  for  low 
frequency  applications  such  as  IR  imaging  systems38. 

The  operation  of  Type  II  photodiodes  under  zero  bias  ensures  that  the  main  noise 
component  is  the  thermal  (Johnson)  noise  and  1/f  noise  is  eliminated.  Our  experimental 
measurements  indicate  that  even  under  a  considerable  reverse  bias,  Type  II  detectors  do 
not  have  a  high  frequency  1/f  noise.  Figure  41  shows  the  frequency  spectrum  of  the 
amplified  output  of  a  detector.  The  device  was  under  a  -0.2  volt  bias  and  illuminated  by 
the  chopped  DR.  radiation  of  a  blackbody.  The  chopper  frequency  was  f0=396Hz, 
blackbody  temperature  and  aperture  diameter  were  Tbb=800  K  and  Dbb=2.54  cm,  and  the 
detector  was  located  d=  15  cm  away  from  the  blackbody  aperture.  Although  the 
measurement  includes  the  noise  of  the  pre-amplifier  (Analog  Device  AD797)  and  the 
FFT  spectrum  analyzer  (Stanford  Research  System  SR  760),  the  knee  of  the  1/f  noise  is 
below  -100Hz.  Under  the  given  parameters  of  the  inset  of  the  Figure  41,  the  signal  to 


noise  ratio  (SNR)  was  more  than  44dB  with  a  bandwidth  of  Af=100  Hz  around 
f0=396  Hz. 


Frequenoy(Hz) 


Figure  41.  Frequency  spectrum  of  the  amplified  output  of  a  detector  illuminated  by  a 
chopped  IR  radiation  of  a  blackbody.  The  detector  was  biased  at  -0.2  V  and  was  located 
d=  1 5  cm  away  from  the  blackbody  aperture. 


III.G.4.  Response  Time 

39 

We  used  a  quantum  cascade  laser  (QCL),  operating  at  room  temperature  ,  as  a 
high-speed  source  of  IR  radiation  at  X=5pm  to  study  the  response  time  of  the  uncooled 
Type  II  devices.  Figure  42  shows  the  current  of  the  QCL  as  well  as  the  output  of  the 
pre-amplifier  versus  time.  The  inset  shows  the  schematic  diagram  of  the  measurement 
setup.  The  laser  threshold  current  marks  the  current  above  which  the  laser  starts  emitting. 
Considering  the  fall  time  of  the  pre-amplifier  (~50  ns)  and  the  fall  time  of  the  output 
signal  of  the  pre-amplifier  (~1 10  ns),  the  detector  response  time  was  about 
(1 102-502)1/2=~100  ns  using  the  sum-of-squares  approach. 
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Figure  42.  The  current  of  the  QCL  laser  (right  axis)  and  the  amplified  output  of  the 
detector  (left  axis)  versus  time.  Inset  shows  the  measurement  setup.  Note  that  the  laser 
emits  when  current  is  above  the  threshold  level. 


III.H.  Device  Performance  and  Comparison  with  State  of  the  Art 

The  performance  of  state-of-the-art  uncooled  HgCdTe  and  microbolometers  were 
detailed  in  V.I.  We  demonstrated  uncooled  detectors  based  on  Type  II  superlattices  with 
a  50%  cutoff  wavelength  of  8  |am.  The  detectivity  of  the  devices  with  zero  bias  and 
without  optical  lenses  or  anti-reflection  coatings  was  1.2xl08  cmHz1/2/W.  The  response 
time  of  the  detectors,  measured  with  a  quantum  cascade  laser  at  room  temperature,  was 
about  100  nsec.  The  following  table  compares  the  performance  of  these  devices  with 
available  HgCdTe  and  microbolometer  detectors. 


Parameters 

Type  II  SL 

HgCdTe 

Microbolometer 

Detectivity 
Xc=8pm;  T=300K 

~108  cmHz1/2/W 

~108  cmHz1/2/W 

~108  cmHz1/2/W 

RoA 

Xc=8p,m;  T=300K 

1.4x1 0‘2  Qcm2 

~lxlO'4Ocm2 

Uniformity 

Very  Good 

Poor 

Very  Good 

Zero  Bias  Operation 

Possible 

Possible 

Not  Possible 

Response  Time 

~10’8  second 

~1CT8  second 

~10'3  second 

Table  II.  The  performance  of  uncooled  Type  II  superlattice  photodiodes  compared  with 
the  uncooled  HgCdTe  and  microbolometer  at  8  pm. 
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